Multi-band Gutzwiller wave functions for general on-site interactions 



J. Biinemann and W. Weber 
f~^ ' Inst. f. Physik, Universitdt Dortmund, D-44221 Dortmund, Germany 

o\ _ 

^ : F. Gebhard 

^ \ Inst. Laue-Langevin, B. P. 156x, F-38042 Grenoble Cedex 9, France 

Q ■ (Accepted for publication in PRB 57 (March 15, 1998)) 

Abstract 
U. 

"^ , We introduce Gutzwiller wave functions for multi-band models with general 

j^ \ on-site Coulomb interactions. As these wave functions employ correlators for 

the exact atomic eigenstates they are exact both in the non-interacting and 
__i - in the atomic limit. We evaluate them in infinite lattice dimensions for all 

P^ ■ interaction strengths without any restrictions on the structure of the Hamil- 

o 



^ 



tonian or the symmetry of the ground state. The results for the ground-state 
energy allow us to derive an effective one-electron Hamiltonian for Landau 
quasi-particles, applicable for finite temperatures and frequencies within the 



K*" ' Fermi-liquid regime. 

o 



As applications for a two-band model we study the Brinkman-Rice metal- 

CN ■ to-insulator transition at half band-filling, and the transition to itinerant fer- 

^_| \ romagnetism for two specific fillings, at and close to a peak in the density of 

f^~. I states of the non-interacting system. Our new results significantly differ from 

^^ ■ those for earlier Gutzwiller wave functions where only density-type interac- 

j^ ' tions were included. When the correct spin symmetries for the two-electron 

states are taken into account, the importance of the Hund's-rule exchange in- 

__i , teraction is even more pronounced and leads to paramagnetic metallic ground 

P^ \ states with large local magnetic moments. Ferromagnetism requires fairly 

^ ' large interaction strengths, and the resulting ferromagnetic state is a strongly 

correlated metal. 

'O ; 71.10.Fd, 71.30.-hh, 75.10.Lp, 75.50.Cc 
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I. INTRODUCTION 

In transition metals and their compounds the electrons of the open d shells participate 
in the itineracy of the valence electrons. At the same time their atomic correlations, as 
described, e.g., by Hund's first and second rule, remain important. The competition between 
the electrons' itinerant and local features results in many interesting phenomena: magnetism 
is most prominent, but also metal-to- insulator transitions, high-to- low-spin changes, orbital 
ordering etc. occur in these materials. 

For insulating compounds it is commonly accepted that their magnetic behavior can be 
described in terms of the spins of localized electrons which are coupled by superexchange 
via the ligands. However, there are conflicting views on the magnetism of the metallic state. 
One line of reasoning which dates back to van Vleck and others assumes rather small charge 
fluctuations around the average (atomic) d^ configuration, and the localized spins remain 
a useful concept also in the metallic state (minimum polarity modeltl). The other school, 
starting from the Hartree-Fock-Stoner theory, treats magnetism within a single-particle 
band theory, i.e., in a completely itinerant limit. In particular, spin-density functional 
theory quite successfully describes the ferromagnetism of the iron group metals, not only 
concerning magnetic moments but also such details as the shapes of complicated multi- 
sheet Fermi surfaces.^ In the spirit of a free-electron theory the spin-density functional 
theory generally assumes a local exchange-correlation potential which is a function of the 
local charge and spin densities. The success of this effective single-particle theory is quite 
surprising since, in the atomic limit, it cannot reproduce the open-shell electronic structure. 

In his seminal work GutzwillerB proposed a variational approach to the problem of itin- 
erant ferromagnetism in the Hubbard model.Q In his many-body trial state atomic config- 
urations with large deviations from the average occupancy could be reduced with respect 
to a Hartree-Fock reference state, depending on the value of the variational parameter. 
Therefore, his approach incorporates both the itinerant and localized aspects of itinerant 
magnetism. Gutzwiller introduced an approximate evaluation of his many-body wave func- 
tion, the so-called Gutzwiller approximation, and concluded from his one-band results that 
itinerant ferromagnetism requires large interaction strengths. Later, Brinkman and Rice 
realizedH that the Gutzwiller approximation contained a transition from a paramagnetic 
metal to a paramagnetic insulator in which all electrons are localized. The Brinkman- 
Rice insulator provides an instructive example for the more general class of Mott-Hubbard 
insulators.Qtl 

During the last decade new analytical techniques were developed which allow for an exact 
evaluation of the single-band Gutzwiller wave function in one dimension,! and in the limit 
of infinite dimensions.™^ For the latter case the results of the Gutzwiller approximation 
were found to become exact. Furthermore, the results of the Kotliar-Ruckenstein Slave- 
Boson mean-field theoryS were re-derived f^r the paramagnetic and the antiferromagnetic 
case.&El In this work we extend the single-band formalism of Refs. |12|,|13|, which provided 
exact results in infinite dimensions for the whole class of Gutzwiller wave functions. For the 
one-band case the Gutzwiller variational approach in infinite dimensions and Slave-Boson 
theoriescil'tJ on mean-field level were shown to be completely equivalent; see Ref. [0, Chap. 3] 
for a recent review. 

The case of multi-band systems poses a more complicated problem. Each lattice site 



represents an atom with an incompletely filled shell. Consequently, the atomic Hamiltonian 
-f^at should include the relevant properties of the electronic structure of isolated atoms or ions. 
This means that it should comply with Hund's first rule, it should reproduce the essential 
features of low-lying multiplet excitations, and it has to incorporate the symmetry of the 
ligand field. The commonly used form for the atomic part H^x, of the multi-band Hubbard 
Hamiltonian includes: (i) an orbital-diagonal density-density interaction of strength U as in 
the single-band case, and an orbital-nondiagonal density-density interaction of strength f/'; 
(ii) two-particle spin-exchange terms of strength J such that the ground state of i^at fulfills 
Hund's first rule, i.e., it exhibits maximum spin for J > 0. Frequently, the spin-exchange 
terms are taken into account only partially. If the (orbital) spin-flip terms are neglected we 
are left with H^^^ which contains only density-density interactions on the atomic sites. 

It should be noted that a while ago atomic Hamiltonians similar to H^^t were studied using 
the "local ansatz" j^t3 a scheme in the spirit of the Gutzwiller method. However, the results 
presented were limited to small interaction strengths. Recently, the Gutzwiller method was 
generalized to treat multi-band Hubbard models with local density-type interactions H^^"^ 
of arbitrary strengths. tZrEJ The essential idea was to evaluate correlators for atomic multi- 
electron configurations made up of spin-orbital product states ('Slater determinants'). This 
was possible since H^^'^^ is diagonal in these configurations. In a first step the Gutzwiller 



approximation was used,EZrt3 later it was shown that these results become exact in the 
limit of infinite dimensions.B2l In Ref. [^ we compared our results with those of previous 



generalizations of the Gutzwiller approximation to the case of degenerate bands. 

However, the frequently used treatment of H^t as discussed above still violates the atomic 
symmetry; for an example, see below. The reason is the incomplete form of the exchange 
interaction. To establish the correct symmetry it is necessary to include all exchange terms 
which result from a (spin-conserving) two-particle interaction, i.e., we will have to consider 
the contributions from up to four different spin-orbits. Then, the proper n-electron atomic 
eigenstates are certain linear combinations of the respective n-electron spin-orbit product 
configurations. As a consequence, the optimum way to generalize the Gutzwiller wave 
function to multi-band systems is the use of correlators for the atomic n-electron eigenstates 
instead of the pure spin-orbit product states. 

In this paper we introduce and evaluate such variational wave functions with atomic 
correlations. Our formulation allows for arbitrary orbital bases, including more than one 
orbital type per representation, i.e., more than one type of s, p, or d orbitals. It also allows 
for an arbitrary number of atomic sites in the unit cell. In the limit of infinite dimension, 
exact results for the ground-state energy are given in terms of an effective single-particle 
Hamiltonian which defines the band structure of correlated electrons. Thus, our theory 
naturally extends to finite temperatures (Fermi-liquid regime) . As an example we apply our 
theory to a two-band model and show that the correct treatment of the atomic correlations 
yields a variety of results which quantitatively and, in some cases, even qualitatively differ 
from those using pure density correlations. 

Our paper is structured as follows. In Sect, y we introduce the multi-band Hubbard 
Hamiltonian with purely on-site interactions. The spectrum of the general atomic Hamilto- 
nian is supposed to be known. Then, in Sect. |T|, we specify the class of Gutzwiller wave 
functions with atomic correlations and give the exact results for the ground-state energy 
in infinite dimensions. For the case of pure density correlations we recover our previous 



expressions. OMEJ In Sect. |^, we discuss the example of two partly filled Cg bands in more 
detail. We study the Brinkman-Rice metal-insulator transition at half band-filling, and 
itinerant ferromagnetism for two generic band-fillings. A summary and conclusions close 
our presentation. Technical details are deferred to the appendix. 

II. HAMILTON OPERATOR 

A. Multi-band Hubbard model 

Our multi-band Hubbard modela is defined by the Hamiltonian 

^ = E ^m'^'CS;-' + E ^^;-t = ^1 + ^-t . (1) 

i,j;cr,(T' i 

Here, cf^ creates an electron with combined spin-orbit index cr = 1, . . . , 2N {N = 5 for 
3d electrons) at the lattice site i of a solid. We do not yet specify a periodic lattice, i.e., the 
sites i may also represent ligand atoms. Therefore, the number of orbitals N also depends 
on the site, N = Ni. To keep the notation transparent we will drop this additional index in 
the following, and we will use the notion "orbital" for spin-orbit states. 

The first term in (|l]). Hi, represents an appropriate single-particle tight-binding Hamil- 
tonian. Crystal field terms are included in the orbital energies tff^ = ei-a-- We may also 

allow non-diagonal crystal field terms tff^ for cr 7^ cr' in case of a sufficiently large orbital 
basis (or a sufficiently low atomic-site symmetry). In this paper we do not include spin-orbit 
coupling and we may therefore assume that the terms t^J are spin-independent quantities 
which only depend on the spatial part of the underlying spin-orbit wave functions. 

In our model (|lD we assume that the electrons interact only locally. Separating the 
density-density interactions we may write the atomic Hamiltonian as two terms, 

^;at= E K'"' n^^^n,;^' + J] JP '"''"''''' ^.,^.A^s^.. ■ (2) 

cr,cr' (cr^cr') ((Ti<cr-z)^(cr3>(T4,) 

The exchange-type second term transfers two electrons from the orbitals 0-3 > 0-4 into the 
orbitals cr^ < cr2. The quantities Ui and Ji represent all possible two-particle (Coulomb) 
interactions compatible with the symmetry at site i. 

B. The atomic problem 

In our variational wave functions we will deal with operators which project onto atomic 
n-electron eigenstates \Ti) with arbitrary number < n < 2N . Altogether we will have 2^^ 
eigenstates. Each of these n-electron eigenstates has the proper symmetry, i.e., they can 
be classified according to irreducible representations of the group defined by the symmetry 
of site i. Part of this classification is according to the total spin quantum number as -f^iat 
commutes with (Sj)^. The problem of classification is treated in detail by many authors;Ellc3 
here we refer to the book of Sugano et al.B^ 

We now suppress the site index, and introduce the following notation for all possible 2^^ 
multi-orbital configurations / and the corresponding multi-electron configuration states \I). 



1. A configuration / is characterized by the electron occupation of the orbitals, 

/ e {0; (1), . . . , (2iV); (1, 2), . . . , (2, 3), . . . (2iV - 1, 2iV); . . . ; (1, . . . , 2N)} . (3) 

Here the symbol in (j^) means that the site is empty. Then, there follow all 2N 
one-electron configurations, all N{2N — 1) two-electron configurations -the sequence 
of numbers in round brackets (cti, cr2, . . .) is irrelevant- and so on up to the 2N elec- 
tron configuration (1, . . . , 2A^). For example, (cti, 0-2) specifies one of the 45 possible 
[Ar] 3d^ configurations of a Ti"''"'' ion in the frozen-core approximation. 

In general, we interpret the indices / in (|^) as sets in the usual sense. For example, 
in the atomic configuration /\/' only those orbitals in / are occupied which are not 
in /'. The complement of J is / = (1, 2, ... , 2N)\I, i.e., in the atomic configuration / 
all orbitals but those in I are occupied. 

2. |/| = n, i.e., the absolute value |/| of a configuration indicates the number n of a 
multi-electron state, 

|0| = 0; |(o-i)| = 1; |(o-i, 0-2)1 = 2; . . . ; 1(1, . . . , 2Ar)| = 2iV . (4) 

3. A multi-electron configuration state ('Slater determinant') is constructed as 

\i _ 
|/) = |cri,cr2,...,cr|j|) = JJc+^|vacuum) (cr„ G /) . (5) 

n=l 

The sequence of electron creation operators in |/) is in ascending order, i.e., cr^ < ctj 
for i < j. When we add an electron to the configuration eigenstate |/) with the help 
of the electron creation operator we obtain the configuration eigenstate |/ U cr) up to 
the fermionic sign function 

fsgn(o-,/) = (/Uo-|c+|/). (6a) 

It gives a minus (plus) sign if it takes an odd (even) number of anticommutations to 
shift the operator c^ to its proper place in the sequence of electron creation operators 
in |J U cr). In general, we define for / fl /' = 



fsgn(/',/) = (JUJ'| 



r 

ciJI). (6b) 



ri = l 



4. The operator which projects onto a specific configuration / is given by 

mi = mij = \I){I\ = Y[h^Y[{l~h^) , (7a) 

where the operators rhi fulfill the local completeness relation 



J]m/ = U. (7b) 

At this point we also define the operators 

710 = 11 ; n/ = JJ^ ficr for |/| > 1 , (8) 

which measure the "gross" occupancy of the atom. The gross occupancy operator 
fii gives a non-zero result when applied to |/') only if I contains electrons in the 
same orbitals as /'. However, I and I' need not be identical because /' could contain 
additional electrons in further orbitals, i.e., only / C /' is required. Each gross (net) 
operator can be written as a sum of net (gross) operators 

hi = 2_. ^v 1 (Sa) 

mi = Y^{~lf\'\hi, . (9b) 

For practical calculations the net operators rhi are more useful than the gross opera- 
tors hi because the former are projection operators onto a given configuration /, i.e., 
rhirhii = 6ij'mi, as can be seen from ([fa|). 

5. For I ^ r we denote J = I Hi', I = JU h, and T = JU h with hnl2 = 0. We want 
to describe the transfer of |/i| = I/2I electrons from the orbitals I2 to the orbitals Ji 
whereas the contents of the other \J\ orbitals remains unchanged. The gross operator 
for the transfer of electrons from I2 to Ji is given by 

With the help of the fermionic sign function (^ the net operator for this process can 
be cast into the form 

mi,p = \I){I'\ = fsgn(J,/i)fsgn(J, J2)[n^- H (^ ' 0]^h,h ■ (H) 

'^e-' o-eJ\(/iU/2) 

All these operators are also defined for |/i| 7^ |/2|. Note the useful relation 

^h,l2^h,h = ^l2,h^Il,h , (12) 

which is easily proven with the help of the Dirac representation of the operators rhiji. 



h 


^x)( n ^ 


(o-nS/i) 


(o-nS/a) 



The configuration eigenstates |/) form a basis of the atomic Hilbert space. The atomic 
Hamiltonian (|]) is Hermitian and only states with the same number of electrons are mixed. 
For the Hamilton matrix (i/at)/,/' = (-^l-f^atl-^') we can find a unitary matrix T such that 

{T)+H,tT = diag(Er) . (13a) 

The atomic eigenstates |r) obey 

\r) = J2Ti,r\I), (13b) 

H,,\r) = Er\T) (13c) 

with 

J]r,rT+,, = V , ^r':/ = ^;,r- (13d) 

r 

Since only configuration eigenstates with the same number of electrons mix, the matrix T 
is block-diagonal with |r| = |/| for each block. 
The atomic Hamiltonian can be written as 

^at = J2 ^r^r , (14a) 

r 

where the projection operators 

rhr = |r)(r| = J2Ti,rmij,T+j, (14b) 

1,1' 

fulfill the local completeness relation 

^rfir = 1 . (14c) 

r 

For r = / = we set T0 g = 1. For |r| = |/| = 1 the atomic Hamiltonian does not fix T/ p, 
and we may choose them to facilitate the evaluation of expectation values for our variational 
wave functions. 

III. GUTZWILLER-CORRELATED WAVE FUNCTIONS 

A. Gutzwiller wave functions with atomic correlations 

Gutzwiller-correlated wave functions are Jastrow-type wave functions, i.e., they are 
written as the many-particle correlator Pq acting on a normalized single-particle product 
state |$o); 

\'^g) = Pg\%)- (15) 



Expectation values with the single-particle state |$o) are denoted by 

0° ^ (6)0 = ($o|0|$o) • (16) 

In general, these expectation values can be calculated easily with the help of Wick's theo- 
rem.LJ In the following we will assume that local Fock terms are absent in \^o), i.e., 

(*o|c+^c..^,|$o) = 5^y{%\c+J,.J<^o) = S^ynl^ . (17) 

This is the case when our orbital basis is sufficiently restricted, i.e., when we use only one set 
of orbitals for each irreducible representation of the group of the site which we consider in 
Hi of eq. ([I|). For cubic symmetry this means that we only consider one type of s and/or p 
and/or d orbitals. In cases of lower symmetry further restrictions are possible; for example, 
in tetragonal site symmetry s-type and d{3z'^ — r'^) orbitals may mix. For the Hamiltonian (|I]) 
we thus choose a basis where the orbitals are not mixed locally, i.e., tff^ = for ct 7^ cr'. 
In the appendix we treat the general case without these restrictions. 

The one-particle state |$n) is usually chosen as the ground state of an effective one- 
particle Hamiltonian Hf. Apart from the simplest cases i^f is not identical with H,- in 
general Hf^ has a lower symmetry than Hi. In these cases the restriction ([1^) may also fail. 

The one-particle wave function contains many configurations which are energetically un- 
favorable with respect to the interacting part of the Hamiltonian. Hence, the correlator Pq 
is chosen to suppress the weight of these configurations to minimize the total energy in (p. 
In the limit of strong correlations the Gutzwiller correlator Pq should project onto atomic 
eigenstates. Therefore, the proper multi-band Gutzwiller wave function with atomic corre- 
lations reads 

\'^G)=pG\^o)=llh,G\%), 
i 

Ag = n ^fr ' = n [1 + (-^^^r - 1) m,;r] = 1 + ^^ (-^^;r - 1) m^.v ■ (18) 

r r r 

The 2^^ variational parameters Aj;r per site are real, positive numbers. For Aj:ro = and 
all other Aj^r 7^ the atomic configuration IFq) at site i is removed from |$o)- 

Expectation values in Gutzwiller-correlated wave functions \^g) will be denoted as 

We will frequently use the expectation values for the atomic eigenstates, mj^r = (^i;r), and 
for the gross 



will frequently use the expectation values for the atomic eigenstates, mj^r 
the gross and net occupancy operators, riij = {fiij) and rrii-j = {rhi-j). 



B. Exact results in infinite dimensions 

Even in the one-band case the evaluation of Gutzwiller-correlated wave functions is 
a difficult many-particle problem; see Ref. |0, Chap. 3.4] for a review. It can be solved 
completely in the limit of infinite dimensions without further approximations. For d ^ 00 

8 



the electron transfer matrix elements between two sites i and ? at a distance li — ?'| = 
X]i=i I'i'i ~ ji\ on a (hyper-) cubic lattice have to be scaled aS 

^3 = \j \-7^A , (20) 



where ti^j is independent of the dimension. In this way the kinetic and the potential energy 
compete with each other for all d. The bandwidth of the electrons stays finite, and we 
may even use the proper d- dimensional density of states Po-.ole) for our calculations. The 
essential simplification in the limit of infinite dimensions lies in the fact that only local 
properties of the wave function are needed for the calculation of single-particle properties; 
see Refs. 0, Chap. 5], ||25|, and ||2^ for recent reviews on the limit of infinite dimensions. 



The class of Gutzwiller-correlated wave functions as specified in eq. (|T8|) can also be 
evaluated exactly in the limit of infinite dimensions. We defer technical details of the 
calculations to the appendix, and merely quote the main result of our work at this point. 
In the limit of infinite dimensions the expectation value of the Hamiltonian (|l]) in terms of 



the Gutzwiller-correlated wave function (pISD is given by 



W= E C''^(dS;<)o + E^^;-^^;- + E^^;rm,r, (21a) 

i7^j;(Ti,(Ti i;<T j;r 



It is seen that the variational ground-state energy can be cast into the form of the expecta- 
tion value of an effective single-particle Hamiltonian with a renormalized electron transfer 
matrix t^-"^ . Due to the off-diagonal terms in the local interactions @) the q factors are 
arranged in a non-diagonal matrix qf.^ which determines the quasi-particle bandwidth and 
the strength of band-mixing in the solid state. As shown in the appendix, the elements of 
the "g matrix can be written as 



^0,(1 _ 77,0,) 2-^ \l ^0^0 



n^,(l -n^,) ^ Y mlm^, 



(T I r,r' 



E fsgn(o-',/')fsgn(o-,/) JmO,,^^,)mO,T+(,^^)r(,,u,.'),rrr^,/'T,,r' , (22) 



I.V 



where we suppressed the site index and used the definition (H) of the fermionic sign function. 
Eqs. ( PID and (P^ ) show that we may replace the original variational parameters Xi-r by 
their physical counterparts, the atomic occupancies rrii-r- They are related by the simple 
equation 

^?;r = !¥ ■ (23) 



ml 



i-T 



Due to the local completeness relation the probability for an empty site is a function of the 
other atomic occupation densities. 



mi;0 = 1 - ^ rrii-r - ^ mi-r ■ (24) 

r(|r|=i) r(|r|>2) 

For the inoment we suppress the site index. As shown in the appendix, the parameters Ap for 
atomic configurations with a single electron (|r| = 1) are the eigenvalues of a (2A^) x (2 A'") 
matrix Z whose entries are given by 

Zcr,^' = -^5^,ct' - yZ ^ Yl fegn(cr',/)fsgn(CT,/)T(/u<T'),r7'w/u<T)— ¥^^ • (25) 

'^ r(|r|>2) ^ i{(t,<t'<^i) ('^''^ ) 

The unphysical case m^ = can safely be ignored. The matrix Z is diagonalized by a 
unitary matrix T', 

(T')+Z(r) = diag(A^) (|r| = 1) (26) 

with T^r = T^-p- These entries in the matrix T remained undetermined at the end of 

Sect, g.' 

Finally, the local densities nj-o- can be calculated from (p3) as 



n. 



l\CT 



^ m^./ , (27a) 



/(o-e/) 



where the configuration probabilities for |/| > 1 follow from (|A18b| ) of the appendix, 

2 



mi-i 



EE,/^^«r„-^r '„«,. (27b) 



Hence, all quantities in (^) are now expressed in terms of the variational parameters mj.p 
and the properties of |$o)- 

The remaining task is the minimization of (i/) in (|21| ) with respect to mj;r and |$o)- 
A conceivable though numerically unsuitable way to achieve this goal is the following. For 
fixed mj;r an input wave function |$™) defines local occupancies ^j.'™. The wave function 
|$q"*) is the ground state of the effective one-particle Hamiltonian 

^efr,in^ E ?/'^'"^CS;^' + E^:^^^;'^ + E^^;rm,r, (28a) 

^ = e,.^ . (28b) 

Note that we have to impose the condition that the orbitals do not mix locally in our effective 
one-particle Hamiltonian, tff = 5o-,cr'ei;<T- The local occupancies of \^q^^) serve as input 
for the next step in the iteration procedure. In this way the optimum |$o''*) for fixed mi-r 
is found recursively. After the global minimization of (H) with respect to the parameters 
mi-r the optimum effective one-particle Hamiltonian H^^' °p* defines a quasi-particle band 
structure which is suitable for a comparison with experiments. Furthermore, it can be used 

10 



to derive the low-temperature thermodynamics. Naturally, the application of H'^^' °p* is 
restricted to the description of the low-energy physics (Fermi-liquid regime) .QE^'y 

Another route to finite temperatures is the following. The variational ground-state energy 
is a function of the occupancies in momentum space. Hence it can be used to derive Fermi 
liquid parametersE3 which give access to the low-energy physics of metallic correlated-electron 
systems .iHaS 

C. Gutzwiller wave functions w^ith pure density correlations 

If we ignore the non-diagonal terms j7;'^1''^2,<t3,<t4 ^^ ^^^ atomic Hamiltonian (H) we may 

set 

Tit = Si,r (29) 

in all the formulae of the previous subsection. Under these conditions we recover the 
Gutzwiller wave functions with pure density correlations. In fact, for this class of Gutzwiller 
wave functions the variational ground-state energy is independent of the non-diagonal 
terms j/'^'^i ''^2,0-3 ,0-4 ^j^jg jg ultimately due to the fact that the correlator does not change 
the orbital occupation and that Fock terms vanish in |<l>o) according to ([T7|) . 

Under the condition (pQ]) the Z matrix in (|25|) is diagonal. We use the fact that 



rii-rr = m. 



i:cr 



+ Yl "^^;^ (^°) 



/(|7|>2,o-e/) 

due to @ so that the eigenvalues of Z can be written as 






mla 



Hence, for consistency with (^) we should have 

n^■,. = 4^ (32) 

for Gutzwiller-correlated wave functions with pure density correlations. This result can be 
derived more directly. As shown in the appendix, in infinite dimensions we have 

{fli-cr) = (-Pj;G^i;o-Pi;G)o = (^i;<T-F^;G)o ■ (33) 

For the second step we used the fact that now the Gutzwiller correlator contains density 
operators only; compare (p^S]). With the help of ([A12b|) the result (^) follows immediately. 
Thus, eqs. ( pUD and (^) allow us to express explicitly the probabilities for a single occupancy 
in terms of the local densities in |$o) and the variational parameters rrii-i for |/| > 2. 

For pure density correlations and for wave functions which obey (0) different local 
configurations are not mixed. Consequently, the *q matrix becomes diagonal, as can be 



shown explicitly from (|2^) with the help of (^). For the diagonal elements we find 



11 



/Qi;c 



<.(1 



71- 



y^ ^mi-inii, 



(/U(t) , 



and the expectation value of the Hamiltonian reduces to 



i^i;'^,'^' 



i;I 



where 



U,,i = (/|i^,at|/) = Y. ^^''^' ■ 



(34a) 



(34b) 



(34c) 



o-.cr'e/ 



For translationally invariant systems the above equations (|34f) were first derived in 
Refs. ||T7|JT8[| using a generahzed Gutzwiller approximation scheme. A concise description 



of this semi-classical counting approach can be found in Ref. ||23]; for a mathematically 
well-defined procedure, see Ref. [0. The wave function used in Refs. |l^Jl^,|20[ was defined 
as 



^=n n *;?'■ 

i I{\I\>2) 

The wave functions |\E'o) and |$o) are related by the transformation&il 

2N 



(35) 

(36a) 
(36b) 



<T=1 



Since Psp contains single-particle operators only, both |\E'o) and |$o) are single-particle wave 
functions. The relation between the parameters Qi-i and Xij is given by 

(37a) 
(37b) 



and 



9ifi — 


^i-fi ) 




Aj;(T 


9i;o- - 


Ai;0 




^ijKfi 


9i;I — 


. lcre/-^*;<T 


n'^ - 


|i"|-i 


9i;I 


.Aa-el^i-cr 



\I\ > 2) 



(|/|>2) 



(37c) 



(37d) 



holds in infinite dimensions. 

Indeoendently, the expressions (^4]) were derived by HasegawaE3 and Fresard and 
Kotliail^ who used a generalization of the Kotliar-Ruckenstein Slave-Boson mean-field ap- 
proachllj introduced by Dorin and Schlottmann.o In Ref. |20] we proved that the results of 
these approximate treatments are variationally controlled in the limit of infinite dimensions; 
see Ref. ||2^ for further comparison with previous variational and Slave-Boson mean-field 
approaches to degenerate-band systems. 
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IV. TWO DEGENERATE BANDS 

The formulae derived in the appendix are completely general and apply for all atomic 
Hamiltonians (H) and for all kinds of symmetry breaking in the one-particle wave func- 
tion |$o)- Depending on the complexity of the problem, the numerical treatment of multi- 
band correlations may become rather involved. It appears to be a good strategy to study the 
two-band case first which provides the simplest example of a correlated multi-band model. 
To keep our expressions for the "q and the Z matrix as simple as possible, we chose a simple 
cubic lattice with one atomic site per cell and two degenerate d{eg) orbitals per atom. This 
model should reflect the situation of nickel to some extent. For example, Ni "'""'", e.g. in 
NiO, exhibits two d{eg) holes in a high spin state, and metallic nickel has approximately one 
d hole per site. 

Alternatively, we could have chosen atoms with two P(x,y) orbitals on a square lattice. 
However, such a model might be less meaningful for the study of ferromagnetic transitions 
since, as we will see below, these strongly depend on the structure of the density of states 
which sensitively depends on the dimension. In addition, our formulae become exact for 
Gutzwiller-correlated wave functions in the limit of infinite dimensions, and 1/d corrections 
are expected to be much smaller in three than in two dimensions. From our experience in 
the one-band casety'Ej'El we conjecture that the differences between d = 3 and d = oo are 
actually marginal. 

A. Atomic Hamiltonian 

We label the orbitals d{3z'^ — r^) as 6 = 1 and (i(x^ — y^) as 6 = 2, and introduce the 
spin index a =1, |. There are four spin-orbitals per atom, leading to 2*^ = 16 multi-electron 
configurations (see table |). Then the atomic Hamiltonian reads 



6 cr,cr' a 

+ '^ 2_^^l,cT^2-u^l-u^2,u + ^C( C-|^_|C-^ |C2^^C2 I + C2^|C2 j^Cj^ ^C-^ | 1 . 
u 

For two orbitals, ifat exhausts all possible two-body interaction terms. 

All sixteen eigenstates and their respective energies are given in table |. The one-electron 
states and, due to the particle-hole symmetry, all three-electron states are seen to be degen- 
erate. The only non-trivial case are the two-electron states. The model of two degenerate 
d{eg) orbitals leads to the following restrictions enforced by symmetry: first, as we can 
use real wave functions for d{eg) orbitals, the relation J = Jq holds; second, the relation 
U — U' = 2J follows from the cubic symmetry.^ To see this we address the six two-electron 
states. There is one spin triplet ^^2 with the energy U' — J. In addition, there are three 
spin singlets: one, with symmetry ^Ai, has the energy U + Jq whereas the other two have 
the energies U' + J and U — Jq, respectively. Cubic symmetry requiresc3 that these two form 
the degenerate doublet ^E. This symmetry requirement can be derived by a transformation 
into the equivalent basis |3y^ — r^) and \z'^ — x^). 
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For Gutzwiller wave functions with pure density correlations the exchange-type inter- 
actions (second hne in (^)) do not contribute to the variational ground-state energy. If 
we ignore these terms in (|38|) all configurations |/) are eigenstates of the resulting H\ 



dens 
at 



In addition to the states with zero, one, three, and four electrons as listed in table |, the 
energies of the two-electron states are grouped into three doublets. As listed in table 0, 
there are (i) the two components of the spin triplet with 15*^1 = 1 at the energy U' — J; 
(ii) the states |t, i) and ||, f) at the energy U\ and (iii) the states |0, ][) and |ti,0) at the 
energy U . 

B. One-particle Hamiltonian and density of states 

We will use an orthogonal tight-binding Hamiltonian with first and second nearest neigh- 
bor hopping matrix elements. Furthermore, we apply the two-center approximation for the 
hopping matrix elements and exclude any spin-flip hopping. Then the matrix elements in 
momentum space between the "iz^ — r^ (6 = 1) and the x^ — y^ (6 = 2) orbitals are given 

byii 

ei(fc) = 4l((l/2) cos k^ + (1/2) cos ky + 2 cos k,) + (3/2)4!5(cos fc, + cos ky) 
+41 cos k^ cos ky + [(1/4)45^ + 3^5^] (cos k^ + cos ky) cos k, 
+3t)^^5(cos kx cos ky + (1/4) cos k^ cos k^ + (1/4) cos ky cos k^) , (39a) 

e2{k) = (3/2)41^^^ (cos k^ + cos ky) + ^^^((1/2) cos k^ + (1/2) cos ky + 2 cos k,) 

+41 cos k, cos ky + [(3/4)tii + tfl + (9/4)tii] (cos k, + cos ky) cos k, , (39b) 

and the band-mixing is given by 

ei2ik) = t2i{k) = (v^/2)[-4i + ^kcosfc. - cos/.,) 

+ [(v^/4)4i - v^4i + (3v^/4)42,](cos A;. - cos/c,) cos A;, . (39c) 



Here we put the cubic lattice constant equal to unity. As in Refs. [17,18| the hopping 



parameters have been chosen according to general experience for transition metal energy 
bands, tii = lev, 4^ = 0.25 eV, and t^^jf : t^P,f : t^f ^ = 1 : (-0.3) : 0.1. This choice 
avoids pathological features in the energy bands such as perfect nesting at half filling. 

The one-particle part of the Hamiltonian (|l|) is easily diagonalized in momentum space 
via the transformation 

Vt-x. = cos 0fc c-l.i,a + sin (pk ct.2,a > (40a) 

Vt,2,a = - sin (pk c^;i,„ + cos 0fc 0+2^^ , (40b) 

with 

tan(20,) = .^f^^^\,. . (40c) 

(^i{k) -e2(fc) 

The dispersion relations for the rj bands become 
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The r] bands are degenerate along the hne (^, ^, ^) in the irreducible Brillouin zone, the total 
bandwidth is W = 6.6 eV. Our one-particle state |<l>o) is chosen as 

l^^o) = n n ^1,-^2, J vacuum) . (42) 

The Fermi surfaces of both rj bands are invariant under the symmetry operations of the 
lattice. 

The condition (0) is fulfilled due to the degeneracy of the Cg orbitals. For the same 
reason the projected orbital densities of states 



V,{e) = iYl cos'(0fc)<5(e - E,{k)) + sm^{<f)k)6{e - E^ik)) , 

k 

V2{e) = -^ ^sin2(0,)5(e - Ei{k)) + cos^{<j)k)S{e - E2{k)) , (43) 



have to be identical, ^i(e) = ^2(e) = ^o(e)/2, and 



ni^ = - I deVoie) (44) 



oo 



is independent of the band index, n^^ = n2„ = n^. 

Since we built in the cubic symmetry into our starting wave function |<l>o) and our 
atomic Hamiltonian (|38|) preserves this symmetry, our selfconsistency cycle will not change 
this property. Therefore, we may set si^o- = S2,a = s^ and ti^o- = h,^ = ta for our variational 
parameters; compare table H for the notation. Note that the number of t electrons and 
I electrons need not be the same, i.e., we still allow for band ferromagnetism. 

For the study of the ferromagnetic transition it is helpful to consider the density of states 
at the Fermi energy, Do(-£'F,(t)- This quantity as a function of the band filling fraction n^ is 
displayed in Fig. |l]. Later, we will study the half-filled case, no- = 0.5, in the context of the 
Brinkman-Rice metal-to-insulator transition, and the fillings rifj = 0.29 and n^ = 0.35 for 
ferromagnetism. The case n^ = 0.29 corresponds to a maximum in the density of states at 
the Fermi energy. There we expect the strongest tendency to ferromagnetism. 

In this work we take the viewpoint of the canonical rather than the grand-canonical 
ensemble. This means that we keep the zero of energy fixed for all band fillings. Then, the 
Fermi energy moves as a function of the band filling and is different for the two spin species 
in the case of ferromagnetism. Of course, this does not change the results because we could 
have kept the Fermi energy the same for both bands and shifted the minority band against 
the majority band to vary the magnetization density. 
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C. Variational ground-state energy 



Now we derive the explicit form of the ground-state energy functional (Um for our exam- 
pie. To this end we first show that the matrix Z in (|25|) is diagonal. From table | we see that 
for |r| > 2 the atomic eigenstates are also configuration eigenstates, T/^r = ^7,r- In this case, 
the factors ^(jU(T),r^wjucr') require JUcr = J=JUcr', i.e., cr = cr'. For |r| = 2 we note 
that J = 7 with 7 7^ <t, <t' has to hold. Now that |r) = ^^/l/2[\al, (T2) ± los, (T4)] according 
to table I we see again that cr = cr' must hold since either (7, cr) = (cr^, cr2) = (7, cr') or 
(7, cr) = (cr^, 0-4) = (7, cr') must be fulfilled ((cti, (T2) fl {as, 0-4) = 0). Since the Z matrix 
is diagonal it follows that the matrix T^ p is the unit matrix. 

The eigenvalues of the Z matrix are A^ = m„/m^^. Then, the relation 



TJIt 



Uli 



is fulfilled for all F, I with IFI 



n„ 



d7 



tp — 7/tj 

= |/|. We then find from (|2|) that 

+ t_, + 2t, + / + - {dA + 2dE + dfi 



(45) 



(46a) 



gives the probabilities for a single occupancy in terms of the multiple occupancies which 
serve as our variational parameters; see table | for the notation. The probability for an 
empty site follows from the completeness relation (|14c|) as 



e = 1 - 2n° - 271° + dY + 4 +dl + dA + 2dE + 4t| + 4t^ + 3/ . 
Along the same lines it can be shown that eq. ( |27b|) reduces to 



,rP"^r 



(46b) 



(47a) 



With the help of (^) and (^5[) it then follows that 



n„ 



n„ 



(47b) 



holds in our model. Similar arguments, as were used to show that the Z matrix is diagonal, 
can be employed to show that the matrix "q is diagonal. Furthermore, the degeneracy of the 
orbitals due to the cubic symmetry requires 



Qc 



Qo 



(47c) 



Hence, the dispersion relations of the rj bands are rescaled by the same factor q such that 
|$o) is unchanged, and our self-consistency cycle terminates after a single iteration. Thus the 
optimum |$o) can be chosen from the start. Nevertheless, we still allow for ferromagnetism 
since E-p^a remained undetermined thus far. 

A straightforward calculation gives the explicit form of the q factors. 



(la 



<{l-nl) 



(vC+yi^) 2 (v^ + 2 



'd, 



+ V^-a 



{-<'){-^) 



.77) 



(48) 
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We denote the kinetic energy of the (1, a) and (2, a) electrons in the uncorrelated state |$o) 
by 



eo-,0 



l^F.a 



deeVo{e) . (49) 

With the help of table |I| we may then cast the minimization problem into the form 

^var, atom ^ ^ ^^-^^^ ^ (^/ _ j)(^TT ^ ^U ^ ^0) 

+2(f/' + J)dE + {U + J)dA + {2U + 4f/' - 2J)(t^ +h + f)- (50) 

This expression must be minimized with respect to the eight variational parameters, rfj , dj , 
(i°, (i°, (i^;, (i^, t|, tj^, and / for a given band-filling rio-. In a paramagnetic situation, ra^ = ra|, 



the number of variational parameters is reduced to five by the relations d 
t„ = t, and Qa- = q. Furthermore, the relations S| = si and ej^^o = ^1,0 hold. 

For Gutzwiller wave functions with density correlations we employ eqs. 
notations of table ||. Now the variational problem reads 



'crcr 
t 



d? = d, 



t; 



and the 



E 



var, dens 



J2q.e^,o + {U'-J){dl^ + di^) 

cr 

+2U'ds + 2t/4 + {2U + 4f/' - 2J)(t| +t^ + f) . 
Here the g factors are given by 

.v/7) 



(51a) 






+ v^ v^+ V^ 



d\ 



{-a){-a) 



(51b) 



In this case our variational parameters are dj , dj , d^, dc, t|, t|, and /. The probabilities 
for an empty site e and a singly occupied site s^ are related to the variational parameters 
by 



s„ = nl- [c?r + t-<x + 2t, + / + 4 + c?s 



jTT 



jii 



e=l-2rVi-2n\ + dY + d\^ + 24 + 24 + 4tt + 4t| + 3/ 



(52a) 
(52b) 



The expression (^T]) is identical to the one used in Refs. [p!7| , p!8| , p2 



D. Brinkman Rice metal— insulator transition at half band-filling 

As a first application of our variational treatment we study the Brinkman-Rice metal- 
to-insulator transition. For a single band and a translationally invariant system we recover 
the original Gutzwiller wave function.H In infinite dimensions this wave function at half 
band-filling is known to describe a continuous transition from the paramagnetic metal to 
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a paramagnetic insulator at U = Ubr above which all electrons are localized (Brinkman- 
Rice insulator) .0 It can be shown, though, that the Brinkman-Rice transition at a finite 
interaction strength is the consequence of the large-d limit, i.e., it is not contained in the 
wave function for any finite dimension.Ell Hence, statements on the metal-insulator transition 
based on our variational description must be taken with care. Even in infinite dimensions 
the Brinkman-Rice transition can be concealed by an (antiferromagnetically) ordered phase; 



see Ref. |]5S| for the one-band case and Ref. |5^ for A^ = 2. It should be clear, though, that 
the onset of long-range order crucially depends on the choice of the matrix elements for the 
electron transfer. In general, there is no perfect nesting between the Fermi surface and the 
Brillouin zone such that antiferromagnetism is not expected to set in for small interaction 
strengths. 

For multi-band systems the Brinkman-Rice transition occurs at integer numbers 1 < 
n < 2N — 1 of electrons per atom. For two bands the transitions for n = 1 or n = 3 are 
continuous like in the one-band case. There, our results do not differ much from those given 
in Ref. [|1^ for Gutzwiller wave functions with pure density correlations. Thus we focus on 
the case n = 2, where, in general, the transition is discontinuous in the bandwidth reduction 
factor q. This means that a jump occurs at the Brinkman-Rice transition from the finite 
value gBR in the metallic phase to g = in the insulating phase. 

For n = 2 the dependence of the five variational parameters dt, ds, d^, t, / as a function 
of the interaction strength U is shown in Fig. |^. Here, the value J = 0.2 U {W = 0.6 U) 
was chosen. For U = (independent electrons) the values of all quantities are equal. As 
U becomes larger the spin-triplet double occupancy dt increasingly dominates the other 
multiple occupancies; in particular, this is true close to the jump at Ubr- All multiple 
occupancies are discontinuous at the Brinkman-Rice transition since, in the insulating case, 
all electrons are frozen into local spin triplets, i.e., we have dt = 1/3 and all other multiple 
occupancies are zero. Note that in the case of pure density correlations the dominance of dt 
is less pronounced; compare Fig. 2 of Ref. |T7[] . 



In Fig. ^ the q values are shown as a function of U for various J/U ratios. The singular 
case J = {U = U') differs from the generic situation both qualitatively and quantitatively. 
The Brinkman-Rice transition is continuous only at this point, and values of J/U as small 
as 10~^ produce finite jumps of a significant size. A (realistic) value of U' = 0.8U (J = O.lf/) 
is enough to reduce the critical interaction strength Ubr for the Brinkman-Rice transition 
by a factor of two; see Fig. |^. Only ioi U' = U {J = 0) all atomic two-electron states are 
degenerate in energy. Thus, near the Brinkman-Rice transition, all double occupancies have 
equal weight, both in the metallic and in the insulating phase. Any finite J value will remove 
the degeneracies and reestablish the generic case. At the singular point of 'zero configuration 
widtha the critical interaction strength can be given by an analytical expression, first derived 
by Lu.il 

As seen from Fig. |^ the critical interaction strength Ubr and the size of the q factor 
strongly depend on the size of the Hund's-rule coupling J/U. In Fig. H we display the 
behavior of ^br as a function of U' /U at the corresponding critical interaction strengths 
Ubr- For the Gutzwiller wave function with atomic correlations (GWatom) a maximum of 
q^^ ~ 0.4 for U' < U occurs near U'/U = 0.9 {J/U = 0.05). A shallow minimum of 
q^^ ~ 0-1 is seen near U'/U ~ 0.14 {J/U = 0.43). On the contrary, the same curve for the 
Gutzwiller wave function with pure density correlations (GWdens) increases monotonically 



as a function of J/U towards g^™*^ :^0.6att/'/f/ = 0. 

In the range < U' /U < 1 we always find q'^^^ > gBR- Moreover, we have Ubr > f^BR*^; 
see Fig. ^ As expected, the metalhc state is stabihzed by the introduction of the full atomic 
correlations. Nevertheless, the two values for the Brinkman-Rice transition are fairly close 
to each other, Ubr ^ f^BR^^- Therefore, it is interesting to plot the value of the q factor for 
the case of atomic correlations at f/ = f/BB^*^. It is very similar to the ^br curve for pure 
density correlations (see Fig. ^. This shows that the q factor sharply -yet continuously- 
drops as a function of U in the region J/br^ ^U < Ubk before it jumps from ^br to zero at 
the Brinkman-Rice transition. 

The Brinkman-Rice transition is discontinuous because the "metallic" and the "insu- 
lating" minimum compete for the global minimum of the variational energy function. In 
contrast to the one-parameter minimization problem of the single-band case (and the two- 
band case for J = 0) the metallic minimum does not smoothly develop into the insulating 
one in the presence of more than one atomic energy scale. The size of the g-factor jump 
measures the difference between the variational parameters in the metallic and in the in- 
sulating phase. Small discontinuities imply that the variational parameters of the metallic 
state at the transition are close to those of the Brinkman-Rice insulator. For large ^br the 
metallic and the insulating minimum are well separated in parameter space. 

In the Brinkman-Rice insulator all sites are in the state with lowest atomic energy. In 
the metal higher atomic states are mixed in. The strength of the mixing depends on the 
energy separation between the atomic levels. For example, the three S = 1 configurations at 
energy U — 3J are separated from the singlet ^E at energy U — Jhj 2 J. Therefore, the 5* = 1 
configurations dominate the metallic state near the Brinkman-Rice transition for large J, 
and gBR decreases with increasing J. When J becomes too large, near U' = {J/U = 0.5), 
the energy of the three-electron states 3U—5J is approaching the energies of the two-electron 
states. Thus, the value of t is enhanced at the expense of the dt parameter. Consequently, 
gBR increases again. In the case of pure density correlations, the two l^zl = configurations 
It, t) and li, i) have energies U — 3J not too much lower than the two configurations It, i) 
and II, t) with energies U — 2 J. This leads to a competition of the respective occupancies 
for all J, and ^bii** i^ ^ fairly smooth function of J. 

In Fig. 1^ we display the paramagnetic {U, U') phase diagram at half band-filling. It is 
seen that the additional atomic correlations (GWatom) stabilize the metallic phase for all 
U > U' {J > 0) compared to the result of the density correlations. The figure also shows the 
gain in the variational energy when we use the Gutzwiller wave functions with full atomic 
correlations instead of pure density correlations. The gain is shown for fixed value U = L^b^^ 
as a function of U'/U. It is quite considerable, of the order of 0.1 eV, for realistic values of 
J/U ^ 0.1. 

For J < {U' > U), the metal is less stable in the presence of full atomic correlations. 
Note that in this parameter range the insulating ground state is different for the two varia- 
tional wave functions: a unique atomic ^Ai state with energy U — \J\ versus two degenerate 
Iti, 0) and |0, ti) states of energy U for pure density correlations. As a consequence, in this 
limit the violation of the atomic symmetry leads to a qualitatively different result for pure 
density correlations. 

Finally, in Fig. |^, we display the size of the local spin {{Si) ) = Si{Si + 1) at half band- 
filling. In the Brinkman-Rice insulator we have Si{Si + 1) = 2 {Si = 1) when J > 0. For 
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J < the local spin drops to zero at Ubr since the singlet ^Ai is the local ground-state 
configuration in the Brinkman-Rice insulator. We again focus on J > 0. For non-interacting 
electrons (f/ = 0) simple statistical arguments apply and the local spin is readily found to 
be (1/2)(1 + l/2)(8/16) + 1(1 + 1)(3/16) = 3/4. 

For J = the local spin increases very slowly with U, up to 1(1 + 1)(3/6) = 1 at Ubr and 
above. Recall that for J = the Brinkman-Rice insulator is highly degenerate. For J > 
the weight of the local triplets becomes more and more important towards the Brinkman- 
Rice transition. This leads to local spin values for U = Ubr as large as Si{Si + 1) = 1.55 
for J/U = 0.1 and even Si{Si + 1) = 1.8 for J/U = 0.45. As seen above in Fig. |, the 
increase in the local spin is most prominent close to the Brinkman-Rice transition which 
again demonstrates the drastic change in the multiple occupancies there. 

E. Itinerant ferromagnetism 



The formulae which we derived in Sect. [IV Q equally apply for the case of ferromagnetism. 



In this subsection we allow for a finite magnetization density M per band in the z direction, 

< M = {rib,^ - rn,,{)/2 < M,,t = n/4 . (53) 

In Fig. 1^ the magnetization M is shown as a function of U for fixed J/U = 0.2 {U' /U = 0.6). 
The critical interaction for the ferromagnetic transition, f/p*°™, is about a factor two larger 
than its value f/p ^ obtained from the Hartree-Fock-Stoner theory. The corresponding values 
Up*^^ always lie somewhat below the values for the Gutzwiller wave function with atomic 
correlations. In general, the relation MHF(f^) > M^ensiU) > Matom(f^) holds, i.e., for all 
interaction strengths the tendency to ferromagnetism is strongest within the Hartree-Fock 
theory and weakest for Gutzwiller wave functions with atomic correlations. Furthermore, 
the slopes of M{U) are much steeper in the Hartree-Fock results than in the presence of 
correlations. 

The properties of the ferromagnetic phase strongly depend on the spectrum of the atomic 
two-electron configurations. To further analyze this point we included the case of Jq = 0, 
which changes only the excited two-electron states. A shift of the curve M{U) results towards 
smaller interaction strengths; for a given magnetization density a smaller interaction strength 
is required as compared to the correct symmetry case J = Jq (see Fig. 0). The effect is more 
pronounced when we go to the Gutzwiller wave function with pure density correlations. In 
this case all exchange terms in (RSf) are neglected. Then, even the ground state is modified 



since the atomic spin triplet with 5*^ = moves up in energy into the range of the atomic 
spin singlets. Again, the magnetization curve shifts to (much) smaller interaction strengths. 
Both results indicate how strongly itinerant ferromagnetism is influenced by the atomic 
n-electron spectra. 

In Fig. 1^ we chose the particle density per band to be n/4 = 0.29 (more precisely: 
n/4 = 0.2941), right at the maximum of the density of state curve; compare Fig. |l|. For 
this case there are finite slopes of the M{U) curves at Up, and a "Stoner criterion" for 
the onset of ferromagnetism applies. In Fig. |^ we chose the particle density per band as 
n/4 = 0.35. As seen from the density of states in Fig. |I], the density of states at the Fermi 
energy Po(-^F,t) + ^o(-Sf,j,) ^^st increases as a function of the magnetization density, and, 
therefore, a discontinuous transition occurs from the paramagnet to the ferromagnet. 
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In the case of pure density correlations a second jump in the M{U) curve is observed, 
which is absent in the other two curves. As was discussed in Ref. |l^, this jump is related 
to another feature of the density of states. In the Hartree-Fock theory this feature is too 
weak to be of any significance in comparison to the interaction energy. When the full atomic 
correlations are taken into account, this first-order jump at a finite magnetization density 
disappears due to the enhanced flexibility of the variational wave function. Nevertheless, in 
this range of a strongly varying magnetization density we find rapid variations of the various 
double occupancies, similar to the behavior near the Brinkman-Rice transition for n/4 = 0.5. 

Another remarkable difference between the Hartree-Fock and the Gutzwiller method lies 
in the approach to ferromagnetic saturation. In the Hartree-Fock theory the magnetization 
saturates at U values about 20% to 40% above the onset of ferromagnetism at U^^ . In 
contrast, in the variational approach saturation is reached at about twice the onset value, 
Usa.t ^ 2t/F- However, even when the minority spin occupancies are zero and {Sf") is con- 
stant, the majority spin occupancies s^ and rfj vary with U since the limit of zero empty 
sites is reached only for t/ — > oo. 

The magnitude of the local spin as a function of U is shown in Fig. p. For U —>■ oo each 
site is either singly occupied with probability 2 — n ot doubly occupied (spin S = 1) with 
probability n-1. Hence, ((5^i)^)oo = (3/4)(2-n)+2(n-l) = 5(n/4)-l/2. For the correlated 
wave functions this limit is reached from above since, for U < oo, charge fluctuations first 
increase the number of spin-one sites at the expense of spin-1/2 sites which turn into empty 
sites. A further decrease of U will also activate the singlet double occupancies and higher 
multiple occupancies. Thus, the local spin eventually reduces below ((Si)^)oo- On the 
contrary, Hartree-Fock theory does not give the proper large-f/ limit for the local spin. 
Instead, the Hartree-Fock limit is given by {{Si)"^)^ = (ra/4)(3 -|- n/2). 

The change of ((S*,)^) at Up is only a minor effect within the correlated-electron approach. 
In particular, this holds for the case of atomic correlations, where about 90% of the local 
spin saturation value is already reached in the paramagnetic state. Again, the Hartree- 
Fock results are completely different. There, the local spin sharply increases as a function 
of the interaction strength since the absence of correlations fixes {{Si)'^)^^{U < Up^) = 

In Fig. ^ we display the J-U phase diagram for both fillings. It shows that Hartree- 
Fock theory always predicts a ferromagnetic instability. In contrast, the correlated-electron 
approach strongly supports the idea that a substantial on-site exchange is required for the 
occurrence of ferromagnetism at realistic interaction strengths. For the case n/4 = 0.29, 
the differences between the phase diagrams for the two correlated-electron wave functions 
are minor. Due to the large density of states at the Fermi energy, the critical interaction 
strengths for the ferromagnetic transition are comparably small, and the densities for the 
double occupancies in both correlated wave functions do not differ much. For the larger 
band filling n/4 = 0.35, i.e., away from the peak in the density of state, the values for Up 
are larger and, in the atomic correlation case, the Gutzwiller wave functions can more easily 
generate local spin triplets while keeping the global paramagnetic phase. 

Finally, in Fig. ITUl we display the energy differences between the paramagnetic and 



ferromagnetic ground states as a function of the interaction strength for J = O.IU. For the 
correlated-electron case this quantity is of the order of the Curie temperature which is in the 
range of 100 K — 1000 K in real materials. On the other hand, the Hartree-Fock theory yields 
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small condensation energies only in the range of f/ ~ 4eV; for larger U, the condensation 
energy is of order U. Including the correlation effects we have relatively small condensation 
energies even for interaction values as large as twice the bandwidth {U ~ 10 eV). 

V. SUMMARY AND CONCLUSIONS 

In this work we constructed Gutzwiller-correlated wave functions with atomic corre- 
lations for general multi-band Hubbard models. We evaluated these many-particle wave 
functions in the limit of infinite space dimensions {d = oo) for the general atomic Hamil- 
tonian for all interaction strengths without any restrictions on the electron transfer matrix 
elements between orbitals on the same or on different lattice sites. Within the metallic 
phase the differences between three and infinite dimensions were found to be small for the 
Gutzwiller wave function for a single band.MH^Ea Therefore, we expect that our results are 
well applicable for the case of physical interest. 

Our variational states consist of a Jastrow-type many-particle correlator which acts on an 
appropriate Hartree-Fock single-particle product wave function. The Gutzwiller correlator 
is chosen to modify the occurrence of atomic multi-electron eigenstates iFj) as compared 
to the uncorrelated (statistical) case. Therefore, our trial states are exact both in the non- 
interacting and in the atomic limit, and they incorporate the essential competition between 
local and itinerant features of interacting multi-band systems. 

The atomic single-particle states of appropriate symmetry (spin-orbits) constitute the 
basis for our one-electron Hamiltonian which describes the motion of the electrons through 
the solid and provides the Hartree-Fock wave function. The atomic multi-electron configu- 
rations |/j) are product states ('Slater determinants') made from the spin-orbit states. If the 
(on-site) electron-electron interaction contains only density-type two-particle interactions, 
the configurations \Ii) will not couple, and the local probabilities nii-j of these n-electron 
configurations (n = |/|) can be used as variational parameters to minimize the ground-state 
energy. In general, however, the states |/j) do not exhibit the correct symmetry of atomic 
n-electron eigenstates |Fj). The correct symmetry can be established only when all exchange- 
type terms of the atomic Hamiltonian are taken into account, i.e., all atomic correlations 
must be included from the beginning. Hence, all configurations |Jj) within the subspace 
|/i| = n are coupled, and the diagonalization of the atomic Hamiltonian by unitary matrices 
Tj./r results in the n-electron atomic eigenstates |Fj). Therefore, the multiple occupancies 
rrii-r for the atomic eigenstates are the appropriate variational parameters in our problem. 

In many cases the elements of the matrices T are given by symmetry alone; in general, 
however, they must be obtained from the diagonalization of the atomic Hamiltonian. 

The exact results in infinite dimensions can be cast into the form of an effective single- 
particle Hamiltonian with reduced electron transfers between the lattice sites. Since the 
atomic eigenstates are non-trivial linear combinations of one-particle product states, the 
hopping reduction factors are arranged in a 2Nx2N matrix q^ with spin-orbit indices cr, cr'. 
These quantities are non-trivial functions of (i) the variational parameters mj^r? (h) the local 
occupancies of the Hartree-Fock wave function, and (iii) the one-electron densities in the 
interacting case. The derivation of the latter requires the diagonalization of a 2N x 2N matrix 
Z. Further complexities occur for the most general case, i.e., for an extended spin-orbit basis 
with more than one orbital-type per representation of the symmetry group of the atomic 

22 



site. Here we have nonzero values for the orbital-nondiagonal parts of the on-site one-particle 
expectation values. At the expense of another unitary transformation this most general case 
is also covered by our formalism. Naturally, this leads to a more complicated structure for 
the matrices ^ and Z and the effective local hybridizations and one-particle densities. 

Like the density-functional theory the variational method is intrinsically limited to the 
description of ground-state properties, e.g., the ground-state energy, compressibility, magne- 
tization, and magnetic susceptibility. Similar to the density-functional theory our variational 
approach naturally extends to finite temperatures and low-frequency excitations since our 
variational ground-state energy corresponds to that of an effective one-particle Hamiltonian. 
The "correlated bands" of this Hamiltonian can be used for a comparison with measured 
dispersion curves and effective masses. Note that our approach is completely general and 
applies to all multi-band systems. Therefore, we hope that it will be fruitful for the de- 
scription of correlated electron systems in the metallic phase. Naturally, any quasi-particle 
approach is limited to the region of the validity of Fermi-liquid theory. 

In this work we presented explicit results for a degenerate two-band model as the simplest 
non-trivial application of our method. We assumed e^-type orbitals on sites of a simple 
cubic lattice. For the single-particle Hamiltonian nearest and next-nearest neighbor transfer 
matrix elements were used which give rise to two bands of width W = 6.6 eV. In our model 
we included all possible two-particle interactions. Yet, there exist only two independent 
interaction parameters U, J since the relation U — U' = 2J holds due to symmetry, and, 
likewise, Jc = J is fulfilled for the charge exchange term. For our simple model system the 
Z matrix is the unit matrix and the hopping reduction matrix is diagonal. 

As a first application we studied the Brinkman-Rice metal-insulator transition at half 
band-filling. Above some finite interaction strength all electrons localize. As for the one-band 
case this localization transition is rather questionable as a scenario for the Mott transition in 
multi-band Mott-Hubbard systems. The lattice sites will not be isolated as in the Brinkman- 
Rice insulator but they will remain coupled via the itinerant exchange. Thus, in the large-L'^ 
limit, we should expect an antiferromagnetic insulating ground state for J > whereas, for 
J < 0, antiferro-type orbital ordering appears to be most likely. 

Apart from the singular point J = the metal-insulator transition is discontinuous 
as manifests itself in finite changes gsR of the bandwidth reduction factor at the transi- 
tion. The results for Gutzwiller wave functions with atomic correlations significantly differ 
from those for pure density correlations. In particular, this applies to the behavior of the 
curves qBK{J /Ubk) for large J which monotonously increases (decreases) as a function of J 
for atomic (pure density) correlations for J /U > 0.05. The Gutzwiller wave function with 
atomic correlations is seen to be more "flexible" than that with pure density correlations in 
the sense that the metallic state can much better adapt itself to the insulator. 

The general aspect of a discontinuous metal-insulator transition could be generic for 
multi-band Hubbard models. In the insulator the atoms are dominantly in a specific n- 
electron ground state which is compatible with Hund's first rule. Other n-electron states (of 
excitation energy J), even more nil electron states (excitation energies U), are separated 
from the ground state by finite gaps. In the metallic phase a macroscopic number of ener- 
getically unfavorable nil electron states is created and, consequently, also a macroscopic 
number of the other n-electron states. It appears to be unlikely that all of these occupation 
densities change continuously at a single critical interaction strength. Instead, the metallic 
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state breaks down discontinuously when the gain in kinetic energy can no longer compen- 
sate the intra-atomic gaps. However, variational statements on the nature of the transition 
between the metal and the Mott-Hubbard insulator must be taken with great care.QlJ 

Nonetheless, for reasonably small J /U and for the case of a general band structure 
without the perfect-nesting property we expect a transition to an antiferromagnetic state 
with strong electronic localization, i.e., with large charge-transfer excitation energies. We 
do not expect transitions to an antiferromagnetic metallic or small-gap insulating state. Yet 
to test this conjecture an antiferromagnetic trial state needs to be investigated. 

As a second application we addressed the issue of itinerant ferromagnetism. For this 
purpose we chose two band fillings, the first one at the maximum of the density of states, 
the second one close to it. Again, we find a large "flexibility" of the Gutzwiller wave function 
with atomic correlations: the paramagnetic state accommodates large local spins, as much 
as 90% of the saturation value, and only a small jump is observed at the ferromagnetic 
transition. Hence, the paramagnetic metallic state near the transition and, moreover, the 
ferromagnetic state are highly correlated. 

In general, the ferromagnetic transition is found to occur at fairly large interaction 
strengths f/p, with values 1 < Up/W < 2. In addition, finite values of the exchange in- 
teraction J are required with Jmin ~ O.IU. These results may change towards smaller values 
of (J, U) if a larger value for the peak density of states is chosen. In any case, our results 
stress the importance of the atomic Hund's-rule exchange for ferromagnetism in multi-band 
models, a view fostered a long time ago by van Vleck.El The ferromagnetic condensation 
energy is an estimate for ksTc, the Curie temperature for iron-group metals. It is found to 
be of right order of magnitude, Tq ~ 500 K, for interaction strengths as large as 10 eV. The 
condensation energy is a smooth function of the interaction strength, i.e., the ground-state 
magnetization does not too sensitively depend on U. 

In contrast, the corresponding Hartree-Fock treatment yields completely different results. 
The ferromagnetic transition is predicted to occur for small values of U, Up^ /W < 1. The 
spin exchange J is fairly unimportant and the magnetization saturates almost immediately 
as a function of U. Finally, apart from a small interval above U^^, the condensation energy 
is grossly overestimated. Thus, itinerant ferromagnetism in interacting multi-band systems 
is a correlated-electron problem that cannot be treated within a weak-coupling approach. 

The inferior results of the Hartree-Fock treatment might be taken as an indication that 
spin-density functional theory is also inadequate for the description of itinerant ferromag- 
netism in iron-group metals because our results suggest strong correlation effects there. On 
the other hand, the success of this effective single-particle theory may point to an inadequacy 
of our multi-band model for the following reason. The results from spin-density functional 
theory indicate that the minority spin bands are broader and, accordingly, the correspond- 
ing wave functions more extended in space than those of the majority bands. This "orbital 
flexibility" makes the minority spin density to dominate in the interstitial regions. Orbital 
flexibility is not included in the present form of our multi-band models. If considered, e.g., 
by extending the orbital basis, the required interaction strength for ferromagnetism may 
be reduced considerably towards a less correlated situation. In principle, our treatment of 
generalized Gutzwiller wave functions allows us to incorporate such basis extensions. Work 
in this direction is in progress, and applications of our general formalism to real systems are 
presently under investigation. 
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APPENDIX: EVALUATION OF EXPECTATION VALUES 

In the appendix we sketch the essential steps for the exact evaluation of Gutzwiller- 
correlated wave functions in infinite dimensions. First, we choose a new basis in which 
local Fock terms are absent. Second, we select appropriate expansion parameters for a 
perturbation theory around the limit of zero interactions. As a third step we set up a 
diagrammatic theory for the calculation of expectation values based on Wick's theorem and 
the linked cluster theorem.c3 The clue to the exact solution in infinite dimensions is the 
selection of the expansion parameters. They are chosen in such a way that all higher order 
diagrams vanish in infinite dimensions, and the trivial order gives the exact result. Further 



technical details can be found in Refs. |T2|,p!3|,pO| . In the rest of the appendix we derive 



explicit results for the local multiple occupancies and the interacting one-particle density 
matrix. 



1. Change of the local basis 

For a general |$o) the non-interacting local one-particle matrix C° with the entries 

C°^,y = ('^o|e+-y4vl<^o) (Al) 

is not diagonal. Therefore, we derive the formalism here which covers the most general case. 
Since our diagrammatic approach for the evaluation of Gutzwiller-correlated wave func- 
tions requires that local Fock terms are absent, we need to perform a local unitary transfor- 
mation, 

(A2a) 
(A2b) 
(A2c) 

It diagonalizes the non-interacting local one-particle density matrix 

F,) d°F, = diagK';f ) . (A3a) 

This is always possible because C° is Hermitian. As seen from ( |A3a|) , local Fock terms 
are absent in the new basis, i.e., the non-interacting local one-particle density matrix H'^ is 
diagonal. 
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yF+ F. , 


= ^<t,<t' , 


7 




1 





Hl,^„, = {<^o\K^\,,\'^o) = S.,.'{%\h:^^h,.J<^o) = ^.,.'<^ . (A3b) 



/i,0 



For a given |$o) the transformation matrix Fi is fixed, and the local occupancies n^.^ in the 

new basis are the eigenvalues of C\. 

From now on we work in the new local basis. We suppress the site index for the rest of 
this subsection. The notation of Sect. |IIB| remains essentially the same but each operator 



c+ (c^) has to be replaced by /i+ {hfj)- In the new basis the configuration eigenstates are 
denoted by 

|w ^ 
|7Y) = J^ h]^^ I vacuum) (cr„ G H) , (A4a) 

n=l 

and the atomic eigenstates |r) and their projection operator mr = |r)(r| are given by 

\T) = J2An,r\n), (A4b) 

H 

mr = ^ An,r'mn,'H'^r,H' ■ (A4c) 

H,H' 

The elements of the unitary matrix A are given by 

(7^|/)=det(F^,,^.) , {jieI,(Tjen). (A5a) 

The inverse relation to ( |A5a|) reads 

Tj,r = J2An,r{m). (A5b) 

n 

Again, A is block-diagonal. Eqs. ( |A5| ) are defined for |r| = |7i| > 2. For F = 7Y = we 
set ^00 = 1. The entries in A for |F| = |?i| = 1 can be chosen at our convenience. We 
will specify them such that an exact evaluation of our variational wave functions becomes 
feasible in infinite dimensions; see below. 

After the change of the basis we obtain simple expressions for expectation values in |$o) 
with the help of Wick's theorem.E3 For example, we have 

V^ I A |2 h,0 

H 



m„ 



n<'°n(i-nr), (A6) 



where we used the fact that Fock terms are absent in the new basis; see eq. 
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2. Choice of the expansion parameter 



In this subsection we suppress the site index. We proceed along the derivation outhned 
in Refs. p^ . |T3| , |20[| . First, we express the square of the Gutzwiller correlator Pq in terms of 



the operators for the configuration eigenstates, 

r 



n,n' 



(A7) 



Next, we demand that 



X^ yn,H'mH,H' = Yl 

n^H' 1-1,1-1' {\H\,\'H'\>2) 



^H,H'n-H,H' 



such that 



^a = i+ E 



^H,H'ft-H,H' 



H,H' {\H\,\H'\>2) 



Here, 



n 



HF 



xHF 



H,H — "'H 



n-i 



n 



~HF 



n. 



n^^ = n''^-n^ 



h,0 



for Ti = 71' and 



n^ 



,HF 
H,H' 



n 



■HF 



nm,H2 {J = nnn';n = juni;n' = juHq 



(A8) 



(A9) 



(AlOa) 
(AlOb) 

(AlOc) 



for H 7^ H'; compare Sects. ^I B| and |A 1| . Note that ($o|^h^|^o) = because we subtracted 
the Hartree terms, and all Fock terms vanish in |$o) due to eq. ( |A3 b[ ) . 

The expansion of Pq in ([A9|) is chosen such that at least four lines meet at every internal 
vertex in our diagrammatic expansion; see below. The number of parameters x-h^-h' in (|A8|) 
is less than the number of parameters yn,n' due to the restriction |7i| = \7i'\ > 2, i.e., we 
essentially require 



X0,0 = 

Xrrrr' =0 



(Alia) 
(Allb) 



Alternatively, as follows from 



these 1 + (2A^)^ local conditions can be formulated as 



($0|i^Gl'^0) = 1 , 



(Al2a) 
(Al2b) 
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The first equation follows immediately because we eliminated all Hartree terms from the 



right-hand-side of (|A9|) . For the other (2A^)^ equations (|A12b ) we analyze the case 



first. The operators n^ on the right-hand-side of ( |A9| ) contain at least two Hartree-Fock 
operators (n^^ ~^^^){^%^ ~^^^) which cannot be eliminated completely by a single operator 
77,^. The term with n^y^, for Ti ^ Ti' vanishes because of the Fock terms which transfer 
electrons from Ti' to Ti. According to eq. (|A3b| ) the expectation value of Fock terms vanishes 



in |$o)- For (t ^ cr' we note that eq. ( [A1U| ) requires Ti' = 7Ucr and Ti = "yUcr' to eliminate all 
possible Fock terms. Nevertheless, this contribution still vanishes because of the remaining 
Hartree-Fock operator h^^ for the orbital 7. 

In Sect. [A6| we shall give the explicit solution of the equations (|A12|) in infinite dimen- 



sions. 



3. Diagrammatic theory and simplifications in infinite dimensions 

Since the variational parameters obey X^r = 1 in the absence of interactions, the pa- 
rameters Xi-T-c^T-c' go to zero for vanishing interactions. Therefore, they are suitable for a 
perturbation expansion in which the order of the expansion is given by the number of fac- 
tors X. When we perform this expansion we may apply Wick's theorem for the resulting 
expectation values since |$o) is a one-particle state. As usual in perturbation theory around 
a single-particle statecil the resulting contributions can be represented diagrammatically. In 
our theory the "internal vertices" represent the factors Xi-;H,'H'- Besides these internal ver- 
tices there are also "external vertices" which come from the site dependence of the operators 
O. For example, there are two external vertices at the sites i and j for O = hf.^h.,^,. The 
non-trivial result in infinite dimensions stem from the Hartree contributions at the external 
vertices; see below. 

To obtain an expansion in powers of x we set Pf.Q = 1 + Pj and write 

i k=l ii,...,ik j=ii 

where ii, . . . ,ik specify internal vertices. Here, the prime on the sum indicates that all lattice 
sites ii, . . . ,ik are different when we apply Wick's theorem. Consequently, the "lines" of our 
diagrammatic theory are given by the one-particle density matrix for the single-particle wave 
function |$o) for i 7^ j, 

P-f' = (1 - 6,,){%\h;^Lj%) ^ (1 - 6,,){h^^L^^,)o . (A14) 

Note that we do not have to distinguish between "hole" and "particle" lines because all sites 
are different when we apply Wick's theorem.ll3'EliH2l 

To make further progress we have to apply the linked cluster theorem.c3 Unfortunately, 
the restriction on the lattice sums prevents its direct application. As shown in Ref. |1T3|J20[] 
this problem can be circumvented by a redefinition of the internal vertices, i.e., Xi-^'H,'H' ~^ 
^i;H,H'- As a result we obtain a standard diagrammatic theory with renormalized vertices 
^i;H,H' and lines given by ( |A14|) . Since the trivial order does not contain any internal vertex, 
it is unaffected by the redefinition of the internal vertices. 
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In our theory we subtracted the Hartree contributions and ruled out local Fock terms 
according to (A3b). For our diagrams this implies that there are no trivial loops at any 



internal vertex. Consequently, there are at least three independent paths from one vertex ii 
to another vertex ^2 in each non-trivial diagram since \7i\ = \7i'\ > 2 requires that at least 
four lines meet at every internal vertex; paths are independent if they do not have a line in 
common. In the limit of infinite dimensionsQ only ii = 12 contributes to a diagram if ii and 
12 are linked by (at least) three independent paths. In our case this implies that the diagram 
simply vanishes since a line linking two identical vertices is zero by the definition (|A14|) .EJ'EJB 



Consequently, the trivial order of our expansion gives the exact result in infinite dimension, 
e.g., only the diagram for {h^^h., ,) with a single line survives the limit d -^ 00. The 
remaining task is the calculation of the trivial (Hartree) terms which stem from the external 
vertices. This will be carried out for the local occupancies and the one-particle density 
matrix in the rest of the appendix. 



4. Local atomic occupancies 

In this subsection we suppress the site index. We need to evaluate (mr). As described 
in the previous subsection, the task is readily solved in infinite dimension, 

mr = (mr) = (PG^r^G)o = A^m^ , (A15) 

where we used the definitions (|14b|) and ([18|) . This proves eq. (^). Furthermore, we may 
use this result to show that the condition ( |A12aD is indeed fulfilled. We have 



(<l>o|Pai$o) = (Y, Ar^r)o = 5^mr = 1 , (A16) 

r r 

because the local completeness relation ( |14c|) holds for the correlated wave function in any 
dimension. 

For later use we write Pq in the form 

Pg=Y >^H,H'rnn,n' , (A17a) 

H,H' 

where we defined 

Ah,h' = Yl ^kAtA^.^, . (Al7b) 

r 

In infinite dimensions we then have 

'mn,n' = {'mn^n') = (^g^h,h'^g)o 

= 2^ ^'Hi,H2'^Hz,HA{'^Hi,H2^H,W^Hi,HA)o (Al8a) 

•Hl,W2,-W3,W4 

K 
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where we used ([T2|) in the last step. In particular, for 7i = 7i' we obtain 



El\ |2 h,0 

K. 



(Al8b) 



With the help of ( |A15|) and the definition (|A17b|) eq. ( p7bD follows. 

Now we solve the {2N)'^ equations ([A12b|) . We multiply both sides of ( [X9| ) with h^h^, 
and take the expectation value with respect to |$o)- With the help of eqs. (|14c|) and ( A4c ) 
eq. ( [A12b| ) becomes 



(^iV)o = 5Z^r 5^ An,rA^^^,{h^h^,mn,n')o ■ 



(A19) 



H,H' 



Apparently we may set Ti = J' U cr' [a' ^ J). With the help of ( [A3b|) we then find that 
Ti' = J U (T and, therefore, cr ^ J . We use the definition of the fermionic sign function (^ 
and ( |A15|) which allow us to simplify the above equation to 



z„y = Y,K,T>^l{K,^ 



(A20a) 



|r|=i 



where the entries of the (2A^) x (2A^) matrix Z are given by 



n 
m. 



to^-.-' - E 5 E f^g^^*^'' ^)fegn(o-, :r)A(^u.'),r^ 



/i.,0 



m, 

|r|>2 ^ J{cr,(T'<^J) 



r,{ju<T) h,o 



m 



(<t,,t') 

(A20b) 



We can safely ignore the unphysical case of m'^^ = 0. All quantities in the matrix Z are 
known as soon as we fix |$o) and our variational parameters mr for |r| > 2. Eq. (|A20a]) 
states that a unitary (2A^) x (2A^)-matrix A' with A'^ p = A^ p diagonalizes the Hermitian 
matrix Z, and that Ap > (|r| = 1) are its eigenvalues. 



{AYZ{A') = diag(A2) . 



(A21) 



Therefore, the conditions ( |A12b|) fix the matrix A^^r for |r| = |cr| = 1, and the expectation 
values for the atomic configurations with a single electron are given by r/ir = Ap?7ip. 



5. q matrix 



As in the previous subsection we have to work in the new basis. Therefore, we start the 
derivation of the "q matrix with a unitary transformation. 



\^j;7iS;7i' ~ Z^ ^ini,<^i ^j;(T{,-f[\^i;cnhj.„/) 



(A22) 



a-i.rr'^ 



In Sect. |A3| we showed that the calculation of the interacting one-particle density matrix 
reduces to 
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(hUK<r'J = ('^0|(A;g/^+.,A;g)(p,;g/^,;^;P,;g)|$0) 



(A23) 



in infinite dimensions. There we also showed that only a single line can join the two external 
vertices i and j. This implies 






0-2,<T- 



Yl K-iAnO^ Yl \/d<ii^t-2,72^. 



j;t2''^2 



(A24a) 
(A24b) 



T2,T2 



which proves the general structure of the variational kinetic energy ( |21b|) , 



m= E ^r'i^.A.^0 



'Ti .Ti 



^.. 






(A25a) 
(A25b) 



CTl ,(7-. ,0*2 ,Cr2 



T2,T2 



Recall that the F matrix is the unit matrix when eq. ([17|) is fulfilled. 

From now on we suppress the site index. To derive the explicit form of the "q matrix we 
use Pq in the form (|A17| ) to write 



PohcrPc — 2^ '^Wi,'H2'^-H3,W4^Wi,-H2^<T^W3,H4 • 

'Hi ^"^.2 1^3 ^'H 4: 

We use the Dirac representation of the operators rhfiyii = \Ti.){H'\ and find 

^l)^2 5'?^3,^4 

= 2_^ ^'Hi,'H2^H3,'Hj^g^i(^,'H3)6n2,H3Ucr'rhni,H4 ? 

''^l,W2,W3,'W4 



(A26) 



(A27) 



where we used the definition of the fermionic sign function (j^). Note that cr ^ 7^3 is now 
required. We see that \Hi\ = IH^I + 1, and our arguments presented in the last subsection 
show that we have T^i = Ti' U cr' {cr' ^ Ti') and 7^4 = Ti' in infinite dimensions. Otherwise, 
local Fock terms would appear in the evaluation of ( |A23D . For cr' ^ Ti' we introduce the 
operator 



m-i 



U K U (1--? 

-yeH'Xa' feWW 



(A28) 



which allows us to write (7^3 = Ti) 



PGhlPG = Yht' H H \n'u.'),^nu.)>^n,H'Hr^icT,Ti)kgn{a',Ti')m^,^^, (A29) 

0-' 'H'{cT'^n')'H{(T<^'H) 
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in infinite dimensions. When we compare this expression with ( |A24a|) we see that 



'^^ = Yl 5Z ^(W'u<T'),(wu<T)A7^,w'fsgn(cr,7^)fsgn(cr',7^')(m^'_^,)o , (A30) 

because we finally singled out the electron creation operator h^, for the contraction according 
to Wick's theorem in (|A24a|) . We use eqs. ( |A15|) and ( |A17b|) and the trivial relation 



V^H',H')^ 



to derive the "g matrix in the form 



hfl 



hfl hfl 



-I hfl \ /i,0/-| h,0\ 

1 - n^, \ n^, (1 - n^, ) 



(A31) 



/ 1 ^^ Imrmr' 



(A32) 



X Y^ isgn{a\nysgn{cr,n)^Jm'^^,^J^,^m'^?A+^^^^J^^A^n'ucT'),^A^,^^,AH,^' ■ 

■H,n' 

When eq. (|T^ holds we recover eq. 



6. Local one-particle density matrix 

Finally, we derive an expression for the interacting local one-particle density matrix, 

C^.n'. = K^'J ■ (A33a) 

Note that the local gross occupancies are the diagonal entries of this matrix, n^ = C-y^-y. We 
express this matrix in the new basis. 



cri,cr'. 



(A33b) 



In infinite dimensions the entries of the interacting local one-particle density matrix H in 
the new basis are readily calculated. 



= 2^ A-Hi,-W2-^W3,W4('^Wi,W2^ai^o-i^W3,-H4)o • 



(A34) 



Hi,H2,-H3,'Ha 



In infinite dimensions we may set Tii = Ti.4 = H', 7i^ = HU (t[, and 7^2 = 7Y U cri with 
cri,(T[ ^ H. We then find 
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^<Ti,o-i = J^ isgn{crj_,n)isgn{a[,n)'^Xn',HuaiXnu,T{,n'm'^° 



■H T-c' 



Y, kgn{(T^,n)kgn{a[, H) ^ W -CT^r,Hu^x^7fu<,r' 5Z ^r^',r^r',w'"^J^° • 



mrmr' 

(A35) 

Therefore, the matrix H is known in terms of the variational parameters mr and the prop- 
erties of the one-particle product state |$o)- If the non-interacting local density matrix C° 
is diagonal, i.e., eq. (|l^ is fulfilled, the matrices C and H are identical, and eq. ( |A35| ) also 
leads to <^. 

The matrix H is Hermitian and can thus be diagonalized with the help of the unitary 
matrix X, 

{x\ ^HX = diag(n^) , (A36) 

where n^ are the eigenvalues of the matrix H. The entries of H thus obey 

(72, 0-2 

0-2, cr^ 

We transform back into the representation with c operators and find for the interacting local 
one-particle density matrix in the original basis 

'^Tiai = Y '^?2,7^ Yl ^<^2,72^-Y^,^^-^^i,<T2-S^^^,^i^'Yi,<Ti^<^i,^; • (A38) 

-)'2,72 <Tl,(T^,<T2,Cr2 

Recall that the local gross occupancies are the diagonal entries of this matrix, n~^ = C~^^~^. 

The result ( |A38| ) allows us to cast the local hybridization term in the variational ground- 
state energy into the form 



E ^^"^^(^^.S() = E ^^"^^(2^2S^)o , (A39a) 

71,71 72,72 

where the effective local hybridizations are given by 

f72,7^ = y^ ^7i,7l y^ p p+ , Y x^ „ ^X+, ,F+F, , . (A39b) 



7i,7i 'Ti,o-i <T2,cr2 



This expression simplifies if we assume that there are no local Fock terms already in 
the basis of the c operators. Then, the F matrix becomes the unit matrix. Let us further 
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demand that orbitals with different crystal-field energies do not mix, i.e., t'^''^ = S~^ye-y. If 

our one-particle product state |$o) respects this symmetry the X matrix becomes the unit 
matrix, and we find 



fl2,-y2 



72 



(A40) 



Thus, we recover ( P8b|) for the effective crystal-field energies. 
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FIGURES 

FIG. 1. Model density of states at the Fermi energy as a function of the band filhng ria- = n/4:. 
The dashed hnes indicate the filhngs used in Sect. |^ The total bandwidth is IF = 6.6 eV. 

FIG. 2. Variational parameters as a function of U for J = 0.2[/ ([/' = 0.6C/) at half band-filling. 

FIG. 3. Bandwidth renormalization factor q at half band- filling as a function of U for various 
values of J ( J = {U -U')/2). 

FIG. 4. Bandwidth renormalization factor at the Brinkman-Rice transition U = U-qk as a 
function of U' /U for the Gutzwiller wave function with atomic correlations (full line) and pure 
density correlations (dotted line) . Also shown is the value of the q factor for GWatom at C/ = U^^ 
(dashed line). 

FIG. 5. Phase diagram and critical interaction strength for the Brinkman-Rice transition in 
Gutzwiller wave functions with atomic (GWatom) and pure density (GWdens) correlations as a 
function of U' /U (left Y-axis). The dashed curve shows the energy gain for atomic correlations 
against pure density correlations at [/ = f/gfj'^ (right Y-axis). 



— * Q 

FIG. 6. Size of the local spin {{Si) ) in the paramagnetic Gutzwiller wave function with atomic 
correlations as a function of the interaction strength and various values of J = ([/ — U')/2 for half 
band-filling. 

FIG. 7. Magnetization density per band as a function of U for J = 0.2U for the Hartree~Fock 
solution (HF), the Gutzwiller wave function with pure density correlations (GWdens)) and the 
Gutzwiller wave function with atomic correlations (GWatom) for (a) n/4 = 0.29 and (b) n/4 = 0.35. 
The dotted line indicates the results for GWatom with Jc = 0. 



FIG. 8. Size of the local spin {{Si) ) as a function of the interaction strength for J = 0.2U 
and band-filling n/4 = 0.35 for the Hartree-Fock theory (HF) and the Gutzwiller wave functions 

(GWdens, GWatom). 

FIG. 9. Phase diagram as a function of U and J for the Hartree-Fock solution (HF) and the 
two Gutzwiller wave functions (GWdens, GWatom) for (a) n/4 = 0.29 and (b) n/4 = 0.35; PM: 
paramagnet, FM: ferromagnet. 

FIG. 10. Condensation energy as a function of U for J = 0.2U for the Hartree-Fock theory (HF) 
and the Gutzwiller wave function (GWatom) for n„ = n/4 = 0.29 (full lines) and n^ = n/4 = 0.35 
(dashed lines). 
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TABLES 

TABLE L Eigenstates with symmetry specifications, spin quantum numbers, energies, and 
notation symbols for the sixteen N = 2 atomic configurations 



# 


Atomic eigenstate F) 


1 


|0,0) 


2 


IT,o) 


3 


|o,T) 


4 


U,o) 


5 


|o,i) 


6 


IT,T) 


7 


(IT,i) + U,T))/V2 


8 


U,i) 


9 


(|T,i)-U,T))/V2 


10 


(in,o)-|o,n))/x/2 


11 


(ITi,o) + |o,Ti))/^/2 


12 


IT,n) 


13 


in,T) 


14 


u,n) 


15 


in,i) 


16 


in,n) 



Symmetry S; 



at 



5: 



at 



energy E-p 



prob. 



ai 



3^2 

'A, 

Eg 
Eg 
Eg 
Eg 

Ai 














e 


/2 


1/2 







S| 


/2 


1/2 







S| 


/2 


-1/2 







«i 


/2 


-1/2 







«i 


1 


1 


U'-J 




4'- 


1 





U'-J 




d? 


1 


-1 


U'-J 




4' 








U' + J 




(Ie 








U-Jc 




(Ie 








U + Jc 




dA 


/2 


1/2 


U + 2U' - 


J 


% 


/2 


1/2 


U + 2U' - 


J 


h 


/2 


-1/2 


U + 2U' - 


J 


ti 


/2 


-1/2 


U + 2U' - 


J 


h 








2U + 4U' - 


2J 


f 



TABLE IL Two-electron spin-orbit states with spin quantum numbers, energies, and notation 
symbols for the case of pure density correlations 



# 



Wave function |/) 



S: 



at 



energy Uj 



prob. 



w 

ds 
ds 
dc 
dr 



6 

7 
8 
9 
10 
11 



IT,T) 
U,i> 
U,T) 
IT,i) 
in,o) 

|o,n) 



1 
-1 








U'-J 
U'-J 

u' 

U' 

u 
u 
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